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Here, we study diffusive spin transport in two dimensions and demonstrate that an intrinsic analog 
to a previously predicted extrinsic spin swapping effect, where the spin polarization and the direction 
of flow are interchanged due to spin-orbit coupling at extrinsic impurities, can be induced by intrinsic 
(Rashba) spin-orbit coupling. The resulting accumulation of intrinsically spin-swapped polarizations 
is shown to be much larger than for the extrinsic effect. Intrinsic spin swapping is particularly strong 
when the system dimensions exceed the spin-orbit precession length and the generated transverse 
spin currents are of the order of the injected primary spin currents. In contrast, spin accumulations 
and spin currents caused by extrinsic spin swapping are proportional to the spin-orbit coupling. We 
present numerical and analytical results for the secondary spin currents and accumulations generated 
by intrinsic spin swapping, and we derive analytic expressions for the induced spin accumulation at 
the edges of a narrow strip, where a long-range propagation of spin polarizations takes place. 



I. INTRODUCTION 

Understanding the spin-orbit interaction is essential to 
the development of spintronics and gives rise to various 
spin transport mechanisms. Effects of the spin-orbit in- 
teraction can be intrinsic or extrinsic. Intrinsic effects are 
caused by the spin-orbit interaction in the band struc- 
ture. Extrinsic contributions arise from spin-orbit cou- 
pling at impurities. The spin Hall effect, where a trans- 
verse spin current is generated via a longitudinal charge 
current, is one of the effects resulting from the spin- 
orbit coupling and has attracted much attention, both 
theoretically-^— and experimentally-^— . In a sample, this 
transverse spin current generates opposite spin accumu- 
lations at the lateral boundaries. 

While the spin Hall effect provides coupling between 
charge and spin, another spin-orbit-induced transport 
mechanism has recently been introduced in which only 
spins couple and which emerges even in the absence of 
charge currents. Primary longitudinal spin currents give 
rise to secondary transverse spin currents due to spin- 
orbit coupling at extrinsic impurities. 16 The generated 
secondary spin currents are proportional to the extrinsic 
spin-orbit coupling strength. The effect has been coined 
'spin swapping' because, in its simplest manifestation, 
it interchanges the spin polarization direction and the 
spin flowJ^ It has been suggested that any mechanism 
inducing a spin Hall effect should also give rise to spin 
swapping. However, it has not yet been clear how the 
intrinsic mechanism could produce this effect. 

In this paper, we demonstrate that an intrinsic (Rash- 
ba spin-orbit-induced) spin swapping effect exists in two- 
dimensional diffusive metals and that it is drastically dif- 
ferent from its extrinsic analog. The main distinction 
between these two effects is that the extrinsic effect is of 
the same order as the spin-orbit coupling strength and is 
thus small, irrespective of the system size. In contrast, 
the intrinsic spin swapping effect is large for system di- 



mensions exceeding the spin-orbit precession length, and 
the secondary spin currents generated by this effect are 
then of the same order as the primary spin currents. If, 
however, the system width is small compared to the spin- 
orbit precession length, the effect is small but leads to a 
long-range propagation of spin polarizations closely re- 
lated to the increase of the D'yakonov-Perel spin relax- 
ation time in narrow strips^. Furthermore, the symme- 
try of intrinsic spin swapping is more complex and richer 
than that of the extrinsic spin swapping effect resulting in 
a non-trivial dependence on the relative orientation of the 
injected spin flow and the spin polarization. We present 
numerical and analytical results for the transverse sec- 
ondary spin currents and accumulations induced by pri- 
mary spin currents in two-dimensional diffusive metals, 
and we compare the intrinsic and extrinsic spin swapping 
effects. 



This paper is organized as follows. We first provide a 
review of the previously discussed extrinsic spin swapping 
effect in Sec. [Hi and we compute the spin accumulations 
and spin currents induced by an injected primary spin 
current in a two-dimensional diffusive metal. In Sec. IIII1 
we discuss the intrinsic spin swapping effect, numerically 
evaluate the spin densities and spin currents generated 
through intrinsic spin swapping, and derive analytical 
results for the resulting spin currents far away from the 
lateral edges of a sample. Next, in Sec. IIV1 we treat the 
case of a narrow strip whose width is small compared to 
the spin-orbit precession length and find analytical ex- 
pressions for the spin accumulations at the lateral edges 
of a sample stemming from the intrinsic spin swapping 
effect. In Sec. |Vj we briefly discuss how the spin swap- 
ping effects could be observed in experiment. Finally, we 
give our conclusions in Sec. IVII 
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II. EXTRINSIC SPIN SWAPPING 

First, we review the extrinsic spin swapping effect 
introduced in Ref. HE and present its features in two- 
dimensional diffusive metals in order to compare it to the 
intrinsic spin swapping effect to be discussed later. The 
Hamiltonian of the system under consideration reads 



U(P) = -^f P + Mmp(p) + V ao (p), 

where p = (x,y) is a two-dimensional coordinate, 



Pi k 



is the elastic impurity scattering potential, and 



Vso(p) = -i7 er x VV imp (p) 



(1) 



(2a) 



(2b) 



is the spin-orbit coupling. p i is the position of the i th 
impurity, A is the area, v(k) is the Fourier transformed 
scattering potential, er = (<j x , a y , a z ) T is a vector of Pauli 
matrices, and 7 is the dimensionless spin-orbit coupling 
strength. 

Considering transport in the diffusive limit, the spin 
diffusion equation reads 



d 2 p h 



(3) 



where /& is the b component of the spin density, b € 
{x, y, z}, and l s [ is the spin-flip length. In order to study 
spin transport, one also needs to define the spin current. 
In the leading approximation, while neglecting spin-orbit 
effects, the spin current is given by the spin diffusion 
current jjQ = —Dd a fb flowing along a and polarized 
along b, where D is the diffusion constant. The spin- 
orbit interaction gives rise to additional terms in the spin 
current. To first order in the spin-orbit coupling strength 
7, when there is no charge current giving rise to the spin 
Hall effect, the spin current is^i 



Jab 



. AO) 

Jab 



+ x 



(jba " Sabj 



(0) 



The term proportional to the swapping constant x relates 
the spin polarization to the direction of flow and results 
in the induction of secondary spin currents, i.e., a 'spin 
swapping' effect For example, a primary spin current 
directed along x will induce transverse spin currents that 
arise as follows, 



if b + x, and 



? -(°) 

Jxb 



;(0) 



Jbx-, 



The first of these transformations swaps the current's 
flow direction and its polarization. In general, this causes 
spin accumulations at the lateral edges of a sample, as 
we shall see shortly. The swapping constant is linear in 
the spin-orbit coupling strength and can be calculated 
explicitly— 



X 



27PF. 



(5) 



where p-p is the momentum at the Fermi level and short- 
ranged scattering potentials are assumed. 

Extrinsic spin swapping arises from the additional 
terms in the spin current Q that are proportional to 
X, whereas the spin diffusion equation ([3]) is unaltered. 
Extrinsic spin swapping therefore affects the boundary 
conditions for an unaltered, conventional spin diffusion 
differential equation. We will see later that the spin dif- 
fusion equation for the intrinsic spin swapping effect is 
altered as well, giving rise to a richer class of phenom- 
ena. 
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(a) Spin density and spin 
current polarized along x 
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(b) Spin density and spin 
current polarized along y 
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FIG. 1. (Color online) The extrinsic spin swapping effect 
in a semi-infinite two-dimensional diffusive metal of width 
L = 4Z s f- Shown are the scaled spin densities and the scaled 
spin currents according to Eqs. ^ and (f?]). A primary spin 



current jxS injected at x — in (a) induces a transverse spin 

Note that 
are 



current j yy through the spin swapping effect in (b) 



the secondary spin accumulation and spin currents in (b) 
linear with respect to the small swapping constant \. If, 
stead, a primary spin current j£y polarized along y is injected, 



(4) fy an( f jxy are illustrated by (a) and the resulting secondary 



spin density and spin current (f x and j yx , respectively) only 
differ from (b) by a sign. 



Jyy Jzz- 



In order to compare the extrinsic spin swapping ef- 
fect with its intrinsic analog to be discussed in the next 
section, we first study the spin polarizations generated 
via the extrinsic spin swapping effect beyond the dis- 
cussion given in Ref. [l6l We consider a semi-infinite 
two-dimensional diffusive metal of width L into which 
a spin current j^J directed along x and carrying spins 
polarized along x is injected at x = 0. We assume that 
the injected current is homogeneous along y at the in- 
jection edge. Further, we assume impenetrable lateral 
sample edges such that no spin current flows through, 
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with j y b(y — ±i/2) = for any spin polarization b. 
The spin-orbit coupling at extrinsic impurities generates 
a transverse spin current j yy on length scales larger than 
the mean free path according to Eq. Q. In turn, this 
gives rise to an accumulation of spins at the lateral edges 
of the sample polarized along y that is anti-symmetric 
in the transverse coordinate y. The spin accumulation 
and spin current are plotted in Fig. [TJ Fig. 1(a) shows 



the polarization along x, and Fig. 1(b) shows the polar- 



ization along y. In the two-dimensional case considered 
here, no transformation into spins polarized along z takes 
place. Note that the extrinsic spin swapping effect and, 
therefore, the resulting secondary spin accumulations and 
spin currents are of the order of the small swapping con- 
stant x- Solving the spin diffusion equation with the 
above-mentioned boundary conditions, the accumulation 



of spins at the lateral edges of a sample can be obtained 
analytically and may be probed experimentally, 



f y (y = ±L/2) = 



2 7^ y 

±- J J^A xKl (xl si ), for L»/ Sf , 

7T D 



± 



^ JxxX -x/l. 

2 D 



for L <C ki, 



where K\ is the modified Bessel function of the second 
kind and first order. This coincides with the numerical 
result illustrated in Fig. [TJ If, instead, a primary spin 
current ji° y polarized along y is injected, the resulting 
secondary spin densities and spin currents (f x and j yx , 
respectively) differ only by a sign according to Eq. (j4j) 
and can also be illustrated as shown in Fig. [TJ 




2. (Color online) The intrinsic spin swapping effect in a two-dimensional diffusive metal of width L = 4l s and length 



16i s . A primary spin current injected at x = in (a) induces an oscillating transverse spin current j yy in (b) through 



FIG. 

L x ~ 

coupling with the z components of the spins in |(c)[ The resulting accumulation of y components of the spins at the sample 
edges in |(b)| is a signature of the intrinsic spin swapping effect. Shown are the spin densities and spin current densities according 
to Eqs. j[8J) and ^ on a relative scale for each plot. Note that all quantities are of the same order of magnitude. 



III. INTRINSIC SPIN SWAPPING 

We now elucidate the nature of the intrinsic spin swap- 
ping effect. The Hamiltonian of a two-dimensional metal 
with intrinsic spin-orbit coupling reads as 

h 2 k 2 

H(k) = —-+cr-h k + v{k), (6) 
Zm* 

where m* is the effective electron mass, k is the electron 
wave vector, and v(k) is the Fourier transformed scatter- 
ing potential. We assume Rashba spin-orbit coupling, 19 

h k = (ak y ,-ak x ,0) T , (7) 

where a defines the spin-orbit coupling strength. When 
a is sufficiently small, such that the spin-orbit precession 



length Z s = (am*) 1 is much larger than the elastic mean 
free path, the spin diffusion equation reads2£ 

d 2 p f x -^f x = U x f z , (8a) 

t s is 

9 2 p fy ^fy = Uyf z , (8b) 

d 2 p f z ~^h = -Ud x f x + d y f y ). (8c) 

The spin current is given b y 20 i 21 

j ab = -Dd a f b + yD(8 ab f z - 5 bz f a ) . (9) 

The diffusion equations ([5]) for the case of intrinsic spin- 
orbit coupling are more difficult to solve analytically than 
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for the extrinsic case because the a;, y, and z spin compo- 
nents are coupled. Therefore, we numerically study the 
spin currents and the accumulations of spins resulting 
from intrinsic spin swapping in a two-dimensional sys- 
tem. Before presenting the numerical results, we discuss 
the simple analytical expressions that can be derived for 
the spin accumulations and spin currents induced by in- 
trinsic spin swapping far away from the lateral edges of 
a sample. The problem can also be treated analytically 
for a narrow strip system whose width is small compared 
to the spin-orbit precession length (see Sec. IIV|) . 

We hrst consider a case analogous to that given for ex- 
trinsic spin swapping. A spin current jix = j X x(x = 0) 
carrying spins polarized along the x direction and di- 
rected along x is injected at x = 0. Again, we assume 
that the injected current is homogeneous along y at the 
injection edge and that the lateral edges of the sample 
are impenetrable, i.e., j y b(y — ±L/2) = for any spin 
polarization b. The situation is, to some extent, similar 
to the extrinsic case depicted in Fig. [TJ However, while 
the swapping effect in this scenario is straightforward for 
the extrinsic case, it is much more complex and rich for 
intrinsic spin swapping. As mentioned before, analytical 
expressions can be found for the spin currents and accu- 



mulations far away from the lateral boundaries, at dis- 
tances much larger than l s . In this region, the influence 
of the boundaries is weak, and the expressions approach 
the limit of a system that is infinite in the y direction. 
We thus find that a transverse spin current j yy flowing 
along the y direction carrying spins polarized along y is 
induced, 



jyy ( x ) _ 
AO) " 

Jxx 

—k r x/l„ 



(V2 - 1) cos (hx/l s ) 



3 + V2 
V7 



sin [kix/lg 



(10) 



where k r n = \J 2\/2 =p 1. This is the intrinsic spin swap- 



ping effect. The induced spin current reaches its maxi- 
mum, \jyy(x max )\/ jxJ w 61%, within one spin-orbit pre- 
cession length from the injection edge at x = 0. The 
injected spin current itself decays away from the spin 
current source at x = 0, 



jxx (x) 
Jxx 



_ -k r x/l s 



cos(kix/l s ) -\ j= sin (kix/lgj 



(11) 




FIG. 3. (Color online) The intrinsic spin swapping effect in a two-dimensional diffusive metal of width L = 4/ s and length 

In turn, this leads to 
an oscillating 



L x — 16l B . A primary spin current injected at x — in (b) induces a transverse spin current j yz in (c) 



In (a 



an accumulation of z spins at the sample edges, which is a signature of the intrinsic spin swapping effect 
spin current polarized along x is only generated close to the lateral edges of the system. Shown are the spin densities and spin 
current densities according to Eqs. (JSJ and Qona relative scale for each plot. Note that all quantities are of the same order 
of magnitude. 



While extrinsic spin swapping in general directly cou- 
ples cc-polarized and j/-polarized spins, in intrinsic spin 
swapping the conversion between x-polarized and y- 
polarized spin currents occurs via spins polarized along 
z as can be seen from Eqs. (J5| and In addition, spin 
currents and spin accumulations oscillate as a function 



of the distance from the injection edge. The situation 
is depicted in Fig. [2] for a sy stem with width L = il s 
and length L x = 1QI S . In Fig. 2(a)| we see that the spin 
current carrying spins polarized along x, which is given 
by Eq. (fTT|) in the bulk, as well as the spin accumula- 
tion decay away from the spin current source at x = 0. 
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The x components of the spins are converted to z com- 



ponents, as shown in Fig. 2(c) which in turn gives rise 



to a swapped transverse spin current j yy , shown in Fig. 
2(b) that is polarized along y. In the bulk, this current 
is given by Eq. (fTU)) . We also see that this swapped spin 
current causes an oscillating spin accumulation at the 
lateral edges, which is a signature of the intrinsic spin 
swapping effect that may be probed experimentally (see 
Sec. IIVI for an explicit expression of this spin swapping 
induced spin accumulation in a narrow strip system). 

Next, we turn to the case in which a homogeneous spin 
current jxy — jxy(x = 0) carrying spins polarized along 
y is injected at x = 0. To analyze this situation, we first 
find an analytic expression for the transverse spin current 
induced through spin swapping far from the lateral edges 
of the system. We find that the primary spin current j xy 
is directly transformed into a transverse spin current, 



.(0) 

Jxy 



-2x/l s 



jxy( x ) 

M ■ 

Jxy 



(12) 



that gives rise to an accumulation of z spins at the lat- 
eral edges of the sample (again, refer to Sec. |IV] for an 
explicit expression for the induced spin accumulation in a 
narrow strip system). In contrast to the case of extrinsic 
spin swapping, Eqs. ((5} and © provide a direct coupling 
between the y and z spins, with the resulting spin current 
having polarization along z (rather than x): again, the 
resulting current is of the same order as the primary spin 
current. It is only near the lateral boundaries that spin 
currents polarized along x are generated as well. This 
spin current leads to an oscillating spin accumulation at 
the sample edges. This situation is depicted in Fig. [3) 

In both scenarios of injected spin currents discussed 
above, intrinsic spin swapping is a much stronger effect 
than extrinsic spin swapping. 



IV. INTRINSIC SPIN SWAPPING IN A 
NARROW STRIP 

In this section, we will consider the special case of 
a strip whose width L is much less than l s . This case 
is interesting because, in such a system, a long-range 
spin swapping effect can be realized, such that the spin- 
swapped accumulation can extend far along the strip, 
over a length much greater than l s . This long-range be- 
havior is closely related to the increase of the D'yakonov- 
Perel spin relaxation time in narrow strips^. 

Due to the small parameter L/l s , the spin-swap ping 
problem can be treated analytically. Following Ref. fla . 
we introduce new spin density variables, 



1 

71 



( ± fx - ify) , "00 = fz 



(13) 



In terms of these variables, Eq. ([5]) can be transformed 
into 

(ide + yJy) V+ (ify-^/a) = 0, (14) 



where is a 3- vector (0i,0o, 0_i) T and J i: i £ {x, y, z}, 
are the corresponding 3x3 angular momentum opera- 
tors for spin 1. Using Eq. (jSJ), the boundary conditions 
can be expressed as 



(idx + j- J v) 1 l'\x=o = I > ( 15a ) 
(i d y - Y J ^\y=±L/2 = 0, (15b) 



where I is determined by the spin current injected at 
x = 0, 



/±i = ^(Ti^+ij^)' Z ° = -^-- ( 15c ) 

(16) 



V2D 

The unitary transformation 

= eiM^-iv/U)^ 

further simplifies Eq. (fl"4"l) to 

2 n2 



id x + T J y (y)) 0-^0 = 0, 



(17) 



where J y (y) 



-iJ x (7t/2-2y/l s )j e iJ x (n/2-2y/l a ) and 



the boundary conditions at the lateral edges of the sys- 
tem then read 

dy4>\v=±L/2 = 0. (18) 

The transformed differential equations (fT7|) and the 
boundary conditions (|18[) are exact equivalent represen- 
tations of the original problem. 

For the case of a narrow strip, L <C l s , one can ex- 
pand J y (y) up to second order in y/l s to obtain J y (y) = 
J z +2(y/l s ) J y — 2(y/l s ) 2 J z and consider the last two terms 
in this expression as a perturbation. Due to Eq. (|18[) . the 
solution of Eq. (|17p can be represented as a Fourier ex- 
pansion in sin ((2n + l)Try/L) and cos (2niry / X) , where 
n is an integer. Further analysis reveals that only a term 
uniform in y is relevant because the other Fourier com- 
ponents decay very quickly along the x direction. The 
equation for cf>, that is, <fi averaged over —L/2 < y < L/2, 
can then be derived from Eq. (jTTJ) asi£ 



(ic^J 2 ) 2 + §(2-jf)0 = O, 



(19) 



where T = 2L 2 /3l 2 . The general solution of this equa- 
tion that converges for x — > oo has the form <p±i = 

A±ie ±2ix/h e -VTx/l B and ^ q = Aoe -V2Tx/l^ CQ _ 

efficients A can be found from the boundary condition 

(HSU). 

If we consider a case analogous to that presented for 
extrinsic spin swapping in Sec. [Ill where a spin current 
jxx carrying spins polarized along x is injected, we find 
I±i = T^fjjxx and 7 = 0. Applying the unitary op- 
erator (|16[) to this boundary condition we obtain in the 
leading approximation 



0±i(x = 0) = ± 



i {0) l 



2TD 



Mx = o) = o. 



(20) 
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From this it follows that 

•(0); 

7 _ ^_Jxxh_ e ±2ix/i B - Vfx/h 

V2TD 



b Q = 0. (21) 



Using Eqs. (|16J) and . we finally obtain the spin den- 
sities 



fx 

fy 

fz 



VTd v ; 

,(0) 



r, jxx Is V -VF: 
VTD k 
•(0); 



(22a) 

x ' u sin(2x/; s ), (22b) 
(22c) 



to first order in y/l s - The accumulation of y spins at the 
lateral edges of the narrow strip caused by the intrinsic 
spin swapping effect reads 

fy (V = ±L/2) = T yf ^ e " ^ sin (2 X /l s ) . 

(23) 

Since V^/ls "C 1, the spin accumulation oscillates and 
slowly decreases along x. 

Considering the second case treated in Sec. IIII| where 
a spin current carrying spins polarized along y is 
injected, a similar calculation yields 

fx = 0, 

•(0); _ 
f _ Jxy h n - VaFa/is 
J y — . — e 



2TD 

•(0), 

r _ o Jxy h V — y / 2fz/i. 

U ~ V2TDls 



(24a) 
(24b) 

(24c) 



Again, the spin densities slowly decay along x but, analo- 
gous to the previous discussion, no oscillation takes place. 



V. EXPERIMENTAL OBSERVATION OF SPIN 
SWAPPING 

In order to observe spin swapping, a primary spin cur- 
rent needs to be injected. This can be achieved in a 
two terminal setup where a spin current is electrically in- 
jected into a two-dimensional diffusive metal from a fer- 
romagnetic electrode. 16 As discussed here, spin swapping 
then gives rise to spin accumulations at the lateral sam- 
ple edges that could be detected experimentally, for ex- 
ample, by optical means 1 ^ or by measuring the interface 
voltage at weak contacts between the lateral boundaries 
and ferromagnet o 12 ' 14 ! 15 . However, in such a setup, an 
electric current is present in the system as well and addi- 
tional spin currents therefore emerge from the coupling 
of charge and spin via the spin Hall effect. In a two- 
dimensional system with extrinsic spin-orbit coupling, 
the spin accumulations resulting from spin swapping at 



the lateral sample edges are polarized in-plane while 
those generated by the electric current via the spin Hall 
effect are polarized out-of-planei 2 -! 3 - This makes it possi- 
ble to experimentally distinguish the two effects. On the 
other hand, in a diffusive system with intrinsic Rashba 
spin-orbit coupling, a uniform electric field gives rise to 
a uniform in-plane spin polarization via the Edelstein 
effect (while it does not produce spin currents). 22 ' 23 In 
contrast, the resulting in-plane accumulation of swapped 
spins generated by a primary spin current with in-plane 
polarization is opposite at the lateral boundaries, as dis- 
cussed above. This difference allows to distinguish the 
intrinsic spin swapping effect and the Edelstein effect in 
experiment. 

Another possibility is the use of a non-local geometry^ 4 - 
where the spin swapping effects could be observed in 
a part of the system where there is no charge current. 
There, an electric current is injected from a ferromag- 
netic electrode on top of a diffusive metal towards a sec- 
ond electrode. A tunnel barrier between the electrodes 
and the metal assures that the current is injected uni- 
formly and it optimizes the polarization of the injected 
electrons. A pure spin current is thus generated in the 
system, propagating in opposite direction of the injected 
charge current and away from the electrodes. This spin 
current will give rise to spin accumulations at the lateral 
sample edges through the spin swapping effect that could 
be detected experimentally. 



VI. CONCLUSION 

In conclusion, we have demonstrated that there is an 
intrinsic analog to the extrinsic spin swapping effect in 
two-dimensional diffusive metals with Rashba spin-orbit 
coupling. We found that the intrinsic effect is drasti- 
cally different because it is large for system dimensions 
exceeding the spin-orbit precession length and gives rise 
to secondary spin currents and accumulations that are 
of the same order of magnitude as the injected primary 
spin currents while leading to a long-range propagation 
of spin polarizations in narrow strip systems. In con- 
trast, the extrinsic spin swapping effect is proportional 
to the spin-orbit coupling strength for any system size 
and is therefore small. Moreover, intrinsic spin swapping 
is more complex and richer than its extrinsic counterpart, 
resulting in a non-trivial dependence on the relative ori- 
entation of the injected spin flow and the spin polariza- 
tion. 

We derived explicit expressions for the transverse spin 
currents in the bulk and numerically computed the re- 
sulting spin accumulations at the lateral boundaries. In 
addition, we derived explicit expressions for the spin ac- 
cumulations in a narrow strip when L <C l s and found 
that the exponential decay of spin polarizations along 
the x direction is greatly reduced in such systems. We 
further gave a brief discussion on how the spin swapping 
effect could be observed in experiment. 
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